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We develop an analytical expression for the self-energy of the infinite-dimensional Hubbard 
model that is correct in a number of different limits. The approach represents a generaliza- 
tion of the iterative perturbation theory to arbitrary fillings. In the weak-coupling regime 
perturbation theory to second order in the interaction U is recovered. The theory is ex- 
act in the atomic limit. The high-energy behavior of the self-energy up to order l/E"^ and 
thereby the first four moments of the spectral density are reproduced correctly. Referring 
to a standard strong-coupling moment method, we analyze the limit U ^ oo. Different 
modifications of the approach are discussed and tested by comparing with the results of an 
exact diagonalization study. 



I. INTRODUCTION 

The theoretical understanding of correlated electron 
systems represents a central problem in condensed- 
matter physics. One of the simplest, but non-trivial 
models that describe interacting fermions on a lattice, 
is the Hubbard model ||, ||. It is nowadays studied 
intensively to gain insight into the fundamental mecha- 
nisms being responsible for itinerant magnetism, metal- 
insulator (Mott) transitions and high-temperature su- 
perconductivity. However, with the exception of the 
one-dimensional case Q , an exact solution is not known, 
and a completely satisfactory understanding of its prop- 
erties has not yet been achieved. 

As was shown by Metzner and VoUhardt ||^ , the Hub- 
bard model is simplified considerably in the limit of high 
spatial dimensions d. However, it still remains a non- 
trivial model, the essential properties of which are com- 
parable to those at low dimensions d = 2,3. An exact 
solution of the d = oo Hubbard model or a reliable ap- 
proximation for the entire range of the model parameters 
will thus provide a proper dynamical mean-field theory 
of the Hubbard model in any dimension Q . 

In the limit d = oo the Hubbard model can be mapped 
onto the single-impurity Anderson model (SIAM) sup- 
plemented by a self-consistency condition [|o[ ^, 
Thereby it becomes possible to make use of various 
methods that are available for impurity problems. 

Using the mapping onto an effective impurity prob- 
lem, essentially exact results for the infinite-dimensional 
Hubbard model can be obtained from quantum Monte 



Carlo (QMC) calculations , [O, and exact diago- 
nalization methods (ED) |]13,]lJ|. However, these ap- 
proaches though being exact suffer from severe limita- 
tions. Within ED calculations one is restricted to a small 
number of orbitals, and thus a smooth density of states 
cannot be obtained. On the other hand, QMC yields 
its results for the discrete Matsubara energies or along 
the imaginary time axis. Therefore, it is difficult to ac- 
cess the low-temperature regime where statistical errors 
become important within the QMC method. Further- 
more, to obtain dynamical quantities such as the den- 
sity of states, the analytical continuation to the real axis 
becomes necessary which cannot be performed reliably 
in all cases. For these reasons the application of ap- 
proximate methods for the infinite-dimensional Hubbard 
model still remains necessary. 

The development of approximate methods should be 
guided by exactly solvable limiting cases and other rig- 
orous analytical results available. Such information im- 
poses a number of strong necessary conditions on any 
approximation. Moreover, comparison with the numer- 
ically exact results of ED or QMC techniques allows for 
a judgment of the quality of the approximation a poste- 
riori. 

For the SIAM weak-coupling approaches are known to 
yield reliable quantitative information jTs], |l6|, ITtI fsl. 



As has been demonstrated by Yosida and Yamada |15 
p^ , perturbation theory in U is quite well behaved for 
the symmetric case when expanding around the Hartree- 
Fock solution. 

Based on these findings, Georges and Kotliar |lO| in- 



troduced a method called iterative perturbation theory 
(IPT) for the d = oo Hubbard model at half-filling 
(n = 1). Within IPT the SIAM is solved by means of 
second-order perturbation theory around the Hartree- 
Fock solution (SOPT-HF) for arbitrary hybridization 
functions followed by the self-consistent mapping onto 
the Hubbard model. IPT leads to convincing results 
as has been proven by comparison with ED and QMC 
studies |o[ Furthermore, a direct treatment 

of zero temperature and real energies is possible. The 
method turns out to be superior to the direct applica- 
tion of SOPT-HF to the Hubbard model. The latter 
does not yield a metal-insulator transition, and Fermi- 
liquid behavior never breaks down p^ . 

The success of IPT may be partly due to the fact 
that at n = 1 the perturbational treatment accidentally 
reproduces the atomic limit. Away from half-filling this 
desirable property is lost, and (the naive extension of) 
IPT is known to give unphysical results. 

In a recent paper p3| Kajueter and Kotliar proposed a 
modification of the IPT scheme introducing an interpo- 
lating self-energy for the SIAM (see also: Refs. ||2j, ^). 
The self-energy exactly reproduces the trivially solvable 
limiting cases of vanishing interaction, U = 0, of the 
completely filled and the empty band (chemical poten- 
tial /i 1^ ±oo), and the atomic limit. For small U it is 
exact up to order J7^. The low-energy behavior {E i— s- 0) 
is fixed by imposing the Friedel sum rule [p6| (equiva- 
lently, the Luttinger theorem ||^, |2|]). At half-fiUing 
the approach reduces to the usual IPT. Finally, the self- 
energy has the correct asymptotic form for high energies 
E ^ oo: The first two coefficients within an expansion 
in 1/E are reproduced exactly. Results for the spec- 
tral density and the integrated spectral weight based on 
the approach have been compared with ED calculations, 
and a rather convincing agreement between both meth- 
ods has been observed indeed . 

The correct high-energy behavior may be quite impor- 
tant, since this is closely connected with the moments of 
the spectral density ||2^. The moments are defined by 

M!i7^ = ^J_ E"^Ad,{E)dE, (1) 

where Ada denotes the one-electron spectral density of 
the localized (d) state in the SIAM. The moments can 
be calculated from the Hamiltonian directly ||3^ . With 
increasing m, however, they include equal-time corre- 
lation functions of higher and higher order. This fact 
limits the number of moments that can be used in prac- 
tice for the determination of the high-energy behavior 
of the self-energy. 



In their recent approach, Kajueter and Kotliar could 
take into account exactly the moments up to m — 2. 
The main purpose of this paper is to demonstrate that 
improvement is still possible. Modifying their approach, 
we will show that a self-energy function can be con- 
structed that respects the m = 3 moment additionally, 
while all other mentioned limiting cases are still recov- 
ered as before. Thereby, higher-order correlation func- 
tions enter the theory. As will be shown, these can be 
expressed without further approximations in terms of 
the one-electron spectral density and can thus be deter- 
mined self-consistently. Our analysis stresses the impor- 
tance of this m = 3 moment, especially what concerns 
spontaneous magnetism. 

The interpolating self-energy for the SIAM is exact 
in the case of small U and in the atomic limit; to 
see whether it can be considered as a reasonable in- 
terpolation between the weak- and the strong-coupling 
regime, we furthermore investigate (analytically) the 
limit [/ I— > oo. Contact is made with a standard strong- 
coupling moment method [the spectral-density approach 
(SDA)] H, H, g. The SDA, within the context of 
the Hubbard model, has proven its usefulness in sev- 
eral previous studies. Thereby, we can provide an ad- 
ditional independent justification for the interpolating 
self-energy. 

The theory is evaluated numerically. Different ver- 
sions of the approach are tested, which are all compatible 
with the mentioned limiting cases. Finally, we compare 
with the results of an exact diagonalization study from 
Ref. ||. 

II. MAPPING ONTO SIAM 

To begin with, we briefly recall the procedure by which 
the Hubbard model can be mapped onto the SIAM and 
introduce some basic notations. 

The Hubbard model reads: 

H = Y1 ~ cLcja + \uJ2 n^an^-a . (2) 

We consider a d-dimensional lattice with hopping be- 
tween nearest neighbors. Provided that the hopping in- 
tegrals are scaled appropriately, T^^) — t = t* /{2\/d) 
{t* = const.), a non-trivial model is obtained in the 
limit rf i-^ oo . The (non-interacting) Bloch-density of 
states for the d = oo simple cubic lattice is a Gaussian 
, a semi-elliptic Bloch-density of states is obtained for 
the Bethe lattice with infinite coordination number [ p5[ . 

The basic quantity to be calculated within the model 
is the one-electron Green function 



Its diagonal elements can be written in the form 



(3) 



J — ( 



E-iz-fi)-j:4E) 



dz. (4) 



For d oo the self-energy of the Hubbard model, 
'E,a-{E), becomes k independent or site-diagonal |2^, 

The Anderson model for a single impurity (SIAM) is 
given by: 



ka a 
+ Uudand-a + ^ Vkd (cl^Cka + 



Cdc 



(5) 

In all practical calculations the hybridization strength 
between the conduction band and the localized d-level, 
Vfcd, enters via the the hybridization function which is 
defined by: 



T/2 

A(£;) = '"^ 



E-ek 



Let us also introduce the impurity Green function: 
From its equation of motion we immediately have: 



(6) 



(7) 



n 



E^{ed- m) - A(£; + /i) - ^do{E) 



where T^daiE) is the d-level self-energy. 

While in the d = oo limit of the Hubbard model all 
spatial degrees of freedom are frozen, the local (tempo- 
ral) fluctuations still constitute a non-trivial problem. 
This, however, is equivalent to the SIAM. The Hubbard 
model and the SIAM can be connected by the following 
self-consistency condition P, [lO| : 



model. In the case of a ferromagnetic phase, the hy- 
bridization function must be spin-dependent. Provided 
that the condition is fulfilled, the self-energy of the 
Hubbard model is identical with the d-level self-energy 
of the impurity problem, Y,„{E) = Yida{E). This also 
implies the corresponding identity between the respec- 
tive Green functions: Ga{E) = Gda{E). 

If one is able to solve the SIAM for arbitrary hy- 
bridization functions lS.„[E)^ the following two-step pro- 
cedure for solving the d = oo Hubbard model is sug- 
gested ||, ^, ^ : Given a hybridization function, we cal- 
culate (by solving the SIAM) the self-energy T,da{E) in 
the first step. In the second step, using Eqs. (|) and (|), 
a new hybridization function is generated. Both steps 
are iterated until self-consistency is achieved. 

In the following we concentrate on the first step which 
represents the actual problem. We intend to derive an 
analytical expression for the self-energy of the SIAM 
that respects various exactly solvable limiting cases and 
other rigorous facts available. 

III. SPECTRAL MOMENTS 

The overall shape of the one-electron spectral density, 

AdAE) = --linGda{E + iO), (10) 
n 

is fixed by their low-order spectral moments to a large 
extent. The definition of the moments is given in Eq. (^. 
A completely equivalent but independent representation 
is easily derived using the Heisenberg equation of motion 
for the time-dependent operators in the definition of the 
spectral density. We obtain: 



M 



(m) 



= ([£™c,.,ctj+) 



(11) 



(8) where CO — [O, Hsiam] - denotes the commutator of an 
operator O with the Hamiltonian, and [•••,•••]+ is the 
anticommutator. The straightforward calculation up to 
TO = 3 yields: 



M 
M. 



(0) 

(1) 

da 
(2) 



= 1 



Cd + U (nd-a) 

+ 2edU(nd-a) + U^{nd-a) 



kd 



AAE + ^i) = E-{€d-fi)~^aiE)~h{GAE))-^ . (9) 

A has to be interpreted as an effective hybridization 
function which provides a coupling of the d level to 
the external bath of conduction electrons that simu- 
lates all temporal degrees of freedom in the Hubbard 



M. 



+ edU'^{nd-a){'2 + {nd-a)) + U^{nd-a) 
+ J2^,'d{^k + 2Id + 2U{nd-A) 

k 

+ U^{nd-a){l- {nd-a))Bd-a . (12) 



Here we have defined ed,k = Cd.fc - Bdc = Bda - M 
and: 



Bda — Ed- 



(?^^^cr)(l - {rida)) 



^Vkdicl^CdailUd-a-l)) ■ 

k 

(13) 

We notice that the m = 3 moment includes higher-order 
correlation functions. 

We can use these explicit results for the moments to 
fix the high-energy behavior of the self-energy. For this 
purpose we consider the following representation of the 
d Green function: 



(14) 



Expanding the denominator in powers of 1/E, we get 

h 

m=0 



m=0 

The coefficients in the 1 /E expansion of the self-energy, 

oo 



m=0 



can be obtained by inserting Eqs. (|T^) and (|T^) and the 
analogous expansion ot the hybridization function into 
Eq. §): 



(0) 

(1) 

da 
(2) 



U (rid-a) 

U^iud^a) (1 - {nd-a)) 
U^iud^a) (1 - {Ud-a}) X 
{^Bd-a + C/(l - {rid-a))) ■ 



(17) 



An approximate expression for the self-energy of the 
SIAM should be consistent with this rigorously derived 
high-energy behavior. 

The Hartree-Fock approximation for the self-energy, 



^'■d^iE)^U(nd-a) 



(18) 



only respects the zeroth-order coefficient in the high- 
energy expansion. The zeroth and the first coefficient 
are reproduced by the self-energy 



^'d}iE) = U{nd-a) + 



U^{nd-a)il-{nd-a)) 
E + n-ed~U(l-{nd-a)) 



which is obtained when applying the Hartree-Fock de- 
coupling scheme not at the first but at the second level 
in the hierarchy of equations of motion for the Green 
function. This is analogous to the "Hubbard-I" approx- 
imation within the context of the Hubbard model. 
The simplest form of a self-energy that implies the cor- 
rect expansion coefficients up to order 1 / E^ is given by: 



4a(^^) = C^(^d- 



U^{nd-a){l- {ud-a)) 

E + ii-Bd-a-U{l-{nd^ 



(19) 



.)) ■ 
(20) 

The higher-order correlation functions included in Bda 
have to be determined self-consistently as well as the 
mean occupation numbers (uda)- Both, E^^^ and S(^\ 
are correct in the atomic limit. Bda reduces to ed in this 
case, and S*^^^ just coincides with the self-energy of the 
atomic limit. 

The self-energy (20) is the result found within the 
spectral-density approach (SDA) ||l[ || ||, ||. Ac- 
tually, the SDA is a standard strong-coupling approach 
to the Hubbard model. Its main idea is to start from 
a two-pole ansatz for the spectral density and to fix all 
parameters in the ansatz such that the moments (up to 
m = 3) are correct. Recent investigations of the d — 2 
Hubbard model |8j point out that the results of the 
SDA are identical to those of the Roth approach |39 
and the Mori-Zwanzig projection technique [^0[ All 
methods yield a two-pole structure for the interacting 
Green function. Neglecting quasi-particle damping (for 
damping effects cf. Ref. ^^), this two-pole structure 
can be made plausible for the strong-coupling regime: 
An analysis of Harris and Lange p3| rigorously shows 
that additional structures merely have a spectral weight 
of the order {t/UY or less. It is remarkable that such a 
simple form (|2^) of the self-energy is able to reproduce 
at least qualitatively the correct dispersion of the Hub- 
bard bands. This has been proven for the d ~ 2 case 
by comparisons ^ H with QMC Q ^ and with ED 
calculations |46) on small square Hubbard arrays. We 
thus believe that the step going from the Hubbard-I so- 
lution, Eq. (|T9|), to a solution, Eq. (pO|), that respects 
the m = 3 moment additionally is quite important, at 
least for [/ i— > co. It is even decisive if spontaneous 
magnetism is considered as it is known that Hubbard's 
original solution predicts magnetic order only under ex- 
treme circumstances g. Within the SDA the spin de- 
pendence of Bcr induces a spin-dependent shift of the 
bands which favors magnetism. Indeed, the results for 
the ferro- and antiferromagnetic Hubbard model seem 
to be qualitatively correct |Q. On the other hand, the 
SDA completely fails to reproduce Fermi-liquid behav- 



ior for small U . At half-filling it predicts an insulator 
for each [/ > and thereby is not able to describe the 
Mott transition as well. Apart from the neglection of 
quasi-particle damping, however, the SDA yields plau- 
sible results for [/ 1-^ oo. 

The conclusion from the preceding discussion should 
be that it may be quite important to account for the 
m — i moment, especially what concerns spontaneous 
magnetism. This will guide our search for an approxi- 
mate self-energy in the case of the SIAM which can be 
identified with the self-energy of the d — oo Hubbard 
model via the self-consistent mapping. 

IV. MODIFICATION OF IPT 

Our approach is a modification of the iterative per- 
turbation theory (IPT) @. Within the usual IPT the 
SIAM is solved by means of perturbation theory up to 
second order in the coupling U. The self-energy is given 
by: 



■"da- - ^ y-d-al -r j . (21) 

The first-order term is the Hartree-Fock self-energy; the 
second-order contribution reads: 



^ ' h^JJJ E-x + v-z 

{fix)fi-y)f{z) + f[-x)J{y)J{-z)) dxdydz . 

(22) 

Here f{x) = l/(exp(/3x) + 1) is the Fermi function, and 
(3 — l/ksT. The Hartree-Fock spectral density 



,(HF 

is obtained from: 



1_ 



^r'(^) = -^Im4r(^ + ^0) (23) 



'^'^ ^ ' E-{ed-r^a)-A,{E + fi)-U{nd-a)' 

(24) 

The parameter /Xo- has been introduced for later pur- 
poses. Within IPT we have Jla — 

Following Kajueter and Kotliar [^, we consider an 
ansatz for the self-energy: 



EdAE) = U{nd-a) + 



(SOC) 



(E) 



(25) 



Oct, ba- and /Ltcr are treated as free parameters which will 
be fixed such that the approximation becomes exact in 



a number of limiting cases. It is assumed that (25) pro- 
vides a reasonable interpolation between the different 
limits. 

In Ref. Kajueter and Kotliar determined the pa- 
rameter to get the correct m — 2 moment of the 
resulting spectral density and the parameter be to get 
the correct result for the atomic limit. Here, in contrast 
Oct and 6ct will be fitted to the m = 2 and to the m = 3 
moments. 

This can be performed straightforwardly: We start by 
expanding the denominator in Eq. (p^) in powers oil/ E 



to obtain the high-energy behavior of S])j^ 



(SOC). 



^(SOC),ps^ _ 1 (m) 



(26) 



rn—1 



where the coefficients are given by: 



\(HF) 



1 - (rid-c 



\(hf; 



(27) 



Here 



(r...)(«^) = - j J{E)AfJ\E)dE (28) 
is a fictive (Hartree-Fock) particle number, and 



E 

k 



(nd.)(HF)(l (HF)) 

Vi.d(cLc,.)(»^)(2(n,_.)(HF) _ 1) (29) 



is the Hartree-Fock value of the higher-order correlation 
functions defined in Eq. ( p^). Comparing with the exact 
coefficients given in Eq. (|17[), we notice that the IPT 
self-energy S^^J"""-* does not have the correct high-energy 
behavior away from half-filling. 

From the equation of motion for the Green function 
((cdo-; c^-o-))^^^^ it can be seen that the hybridization- 
induced correlation functions in the definition of B^^^ 

(1(7 

can be expressed in terms of the localized HF Green 
function: 



1 

{E + iO + fi) g'jI^^^ {E + iO)dE. (30) 



The high-energy behavior of the interpolating self- 
energy can be derived from the expansion (|2^) and from 
(|25|). Comparing with the exact coefficients of the \/E 
expansion in Eq. (fl^) again, we have to choose 



{nd-a){l - {rid-a)) 

(nrf_,)(HF)(l_(„,_^; 



(HF) 



and 



Bd-. -II- B^^l^ + + U{1 - 2{nd-a)) 
C/2("d-.)(HF)(l-(n,_,)(HF)) 



6.= 



(31) 



(32) 



to ensure the correct high-energy behavior of the self- 
energy SrfCT and thereby the correct moments of the re- 
sulting spectral density up to to — 3. The result ( |3l| ) 
and (|3^) reduces to the approach of Ref. if Bda and 
B^^'aie replaced by ed- 
it is easily verified that our approach is correct in 
the atomic limit. Setting Vkd = 0, the Hartree-Fock 
spectral density ( p3| ) reduces to a (5 function which allows 
to calculate the second-order contribution ( p2| ) and thus 
the self-energy ( p5| ) immediately. It turns out to coincide 
with the self-energy of the atomic limit (^9|). 

Next, we have to check the weak-coupling limit: Pro- 
vided that the parameter Jia is chosen such that ]la ^ 
^0 = lAu=o for [/ 1-^ (see below), we have {uda) ^ 
("-dcr)^^^^ for U ^ Q. Furthermore, since Bda ^ -B^ct^^ 
as J7 0, we have h^, ^ \/U . Therefore, expanding the 
self-energy (|5|) in powers of J7, we see that it is correct 
up to order indeed. In particular, this implies that 
all Fermi-liquid properties as described in Ref. |28) will 
be recovered for small U at least. 

Finally, the parameter Ji„ has to be fixed. The most 
natural choice is the following: 

Mcr = A* ■ (33) 
Another possibility is due to an approach of Martin- 
Rodero et al. ^ where is determined from the 
condition: 



[ride 



)("^^ = M . (34) 

In Ref. [ p3[ K ajueter and Kotliar imposed the Friedel 
sum rule |26[| as a condition to fix /Ig-. Via the self- 
consistency condition (^ this is equivalent to the Lut- 
tinger theorem p7[ | which in the case of the d — oo 
Hubbard model reads lEsI: 



(35) 



Let us briefly discuss the implications of the differ- 
ent choices. First we notice that in all cases we have 



Jj-a Mo for J7 I— > as it must be to ensure the correct 
weak-coupling behavior of the self-energy. Furthermore, 
the validity of the approach within all other limiting 
cases that have been considered is not affected by the 
condition chosen. 

Inspecting the original derivation in Ref. |2^, we re- 
call that the validity of the Luttinger theorem depends 
on a number of presuppositions. For example, the theo- 
rem holds if perturbation theory applies. Obviously, for 
small U all conditions (|3|)-(^) yield a theory that is 
compatible with the Luttinger theorem up to order 
at least. Here, another supposition is more important, 
namely lmI]cr{E) = at = 0. In particular, this 
implies T = [ pO[ . Therefore, applying the third condi- 
tion ( ^5| ) does not allow to consider finite temperatures. 
On the other hand, the conditions (33) and (Q) do not 
suffer from this difficulty. 

The second condition (|3^) implies (for a constant hy- 
bridization function Acr) that /lo- = Mo + U {nd-a)- This 
exactly compensates the energetic shift of the Hartree- 
Fock Green function (H) by U{nd^a)- cfj^ thereby 
becomes independent of f/, and sj^^'^'^'' reduces to the 
second-order contribution from a weak-coupling theory 
in which the free [U — Q) instead of the HF d Green func- 
tion is used in the calculation of the proper irreducible 
self-energy. This, however, has not to be confused with 
the plain (or conventional) weak-coupling theory, where 
the free chemical potential /io would have to be re- 
placed by M additionally and which artificially breaks 
particle-hole symmetry |Q. The second-order contri- 
bution within SOPT-HF is recovered only when using 
the first condition (33). Arguments in favor or against a 
particular weak-coupling approach have previously been 
developed by demanding the correct high-energy behav- 
ior (up to order \/E) jl^, |l^. Such reasoning, however, 
is not meaningful in the present context since the correct 
high-energy behavior (up to order l/E"^) is reproduced 
in each case. 

Let us mention that choosing Jl^r — fJ- (first condition) 
introduces a slight complication: Because of the shift of 
^dCT^^ by the HF contribution U{nd-a), it may happen 
that (uda)^'^^^ = (or (nd^)(HF) = i) for strong U and 
T = 0, which means that the parameters and feo- 
are no longer well defined. However, within the limit 
(rid-cr)'-^^-' 1-^ (or 1-^ 1) we have: 



-.(SOC) 



{rid- 



(hf; 



(1 



{rid- 



\(HF) 



E - X + y 



dx dy dz . (36) 



Therefore, although and ba- diverge, the interpolating 
self-energy remains finite. 

So far it can be concluded that there are no differ- 
ences between the three possibilities considered that are 
crucial. However, we notice that the third condition im- 
plies a restriction of the theory to zero temperature. We 
defer further discussion to Sees. VII and VIII. 

V. SELF-CONSISTENT DETERMINATION 
OF HIGHER CORRELATION FUNCTIONS 



E 

kp 



VkdVpdicl^Cpa) = --7-Im / 

J-oo+iO 



fiE)x 



A^{E + ti) [A,{E + n)GdaiE) + h] dE. (41) 

Applying once more the general spectral theorem, per- 
forming a Fourier transformation to time representation 
and using the Heisenberg equation of motion, the third 
term can be written as: 



The definition of the parameter involves the band- 
filling (nda) and the higher-order correlation functions 
included in Bda which are still unknown. A satisfactory 
theory cannot be constructed unless it is possible to de- 
termine these correlation functions without further ap- 
proximations. No problems are introduced for the band- 
filling which may be expressed in terms of the spectral 
density: 



(rida) 



f{E)AaAE)dE. 



(37) 



In the following we demonstrate that also Bda can 
be reduced to the spectral density or the Green func- 
tion, respectively. According to its definition, Bda 
partly consists of a sum of one-particle correlation func- 
tions. Applying the general spectral theorem and 
exploiting the equation of motion for the Green func- 
tion ((cd^; c[-o-))' '^^ h.a.ve: 



^Vkdicl^Cda 



1 



oo+zO 



-Im / fiE) X 

— oo+iO 

A,{E + ^i) Gda{E) dE . (38) 



Now we are left with the higher-order correlation func- 
tions only. Using the commutator 

[Cda ,H]^ =edCda + Ucdand-a + ^ Vpd Cpa , (39) 

P 

the remaining terms in Bd,7 can be written in the form: 



Vkd{cl„Cd<jnd-a) = -^Y'^kd {ci^Cda) 
k k 

]j J2 VkdVpdiclcp,) + ^^]-> -(40) 



kp 



The first term on the right-hand side has just been 
treated above. Using the spectral theorem and the equa- 
tion of motion, one gets for the second one: 



Y.Vkd{cl[cda,H]^) ^J2^kd (-^) Im 



OC pOC 



— oo ^ — oo 



{{cdAt);cl{md{t^t') dE. 



(42) 



Integration by part and back-transformation to energy 
representation finally yields: 



1 



OO + lO 



Y,Vkd{cl[cda,HU = -— Im / f{E) 

E A^iE + ^i) GdaiE) dE . (43) 
Combining all results, we obtain: 



{ndcr){l - (nda)) Bda ('^da)(l " (nda)) ^d 



1 nOO + lV 

-rim fiE)A,{E + n)x 

J -QO+iO 



-JldAE)-!] GdAE)dE, (44) 



where once more we exploited the equation of motion 
for Gda- This completes the theory since (ua-) and Bda 
can be determined self-consistently from the d Green 
function. 

VI. STRONG-COUPLING LIMIT 

So far we have shown that the appropriate choice of 



the parameters and b^ in the ansatz (25) yields a 
self-energy that is exact in a number of limiting cases, 
namely trivially for U — 0, for (nda) — and (nda) ~ 1 
and furthermore in the atomic limit and for small U up 
to order [/^. The high-energy expansion yields the cor- 
rect coefficients up to order (l/E)^. Thereby the first 
moments of the spectral density up to m = 3 are repro- 
duced exactly which is one additional moment compared 
with the approach of Kajueter and Kotliar pSf. 



It is easily shown that the general particle-hole sym- 
metry requires Bda to be zero at half-filling in the para- 
magnetic phase (and for a symmetric Bloch-density of 
states). Since in this case b'"^^ = and fi = Jla- = U/2, 
it follows that ba — 0. Consequently, the theory reduces 
to the conventional IPT for (rida) = {nd-a) = 0.5. 

We now turn over to the strong-coupling case U t—^ oo 
in order to check whether the theory can be regarded 
as a reasonable interpolation between the weak- and 
the strong-coupling regime. In the following we exclu- 
sively focus on the SIAM; any approach that uses the 
self-consistent mapping of the Hubbard model onto the 
SIAM can only be as reliable as the approximation em- 
ployed for the solution of the latter. Hence, for the 
present purposes the hybridization function A can be 
assumed to be independent of U. 

Let us consider the first condition ( p3| ) for determin- 
ing the Active chemical potential, i. e. Jl^ = /i. As 
mentioned above, for small U all Fermi-liquid proper- 
ties are recovered. For large U the spectral density Ada 
is expected to consist roughly of two dominant features 
separated by U at about and + U. Assuming 
{n-da) = (nd-a) < 0.5 for simplicity, the chemical po- 
tential will be located within the lower part: fi ed- 
Looking at the HF Green function (|2^), we notice that 
(ud-cr)^-^^^ = for T = and for U larger than a certain 
critical value Uci, and thus the relation ( ^6| ) applies. In 
particular, at the Fermi edge E = 0: 



.(SOC) 



~'d<7 



(0) 



^^"^ (nrf_,)(HF)(i-(n^_^)(HF)) 
5// i-)ATZ' iy)ATZ^ (y - dy . (45) 

A simple analysis shows this two-fold convolution in- 
tegral to be non-zero for interactions U smaller than 
another critical value Uc2 (> C/ci)- Via Eq. ( p5| ) this 
implies a non-zero value of the imaginary part of the 
interpolating self-energy a,t E = 0. For U > Uc2 we 
have lmYidcr{E) = in the vicinity oi E = {). Therefore, 
non- Fermi- liquid behavior is implied for all U > Uci- 

On the contrary, at T = the second condition (jsj) 
yields liaYida{E) ^ E^ ior E ^ Q irrespective of the 
value for U . 

Some more statements can be made if J7 i—s- oo. For 
this purpose we first determine the total weight and the 
center of gravity of the imaginary part of the second- 
order contribution y}-^^'^^ . Let us consider the integral 



f 

J — c 



( - i ) Iml]|,f {E + iO) dE , (46) 



for m = 0, 1. Since S^^*^*^-* is analytical in the upper 
half of the complex i?-plane, the contour along the real 
axis can be deformed into a contour C where each point 
of C has a distance from that is larger than R and 
i? I— > oo. On the contour C one is allowed to replace 
T!f^'^^ by its asymptotic form (p6|). After that C can 
be re-deformed into an integration along the real axis 
and the evaluation of the integral becomes trivial. The 
total weight and the center of gravity turn out to be 
D^^^ and D^da l^^d}^ respectively, where -D^™' are the 



coefficients from the expansion (|2q). Analogously, the 
total weight and the center of gravity of —lT[vYida{E)/TT 
are given in terms of its 1/E expansion coefficients as 

At energies with \E — D^^^ / D^^^\ i— > oo, the imaginary 
part of the second-order contribution vanishes, and the 
real part approaches: 



ReE 



(SOC) 
da 



(E) 



D 



(1) 



E-D^''/D^'> 

da I da 



(47) 



Now let us assume that there is a zero E = E^'^ of the 
denominator in Eq. (EE 



(48) 



l-5.4f«'(4°)) = 0, 
at a sufficiently high energy, i. e. we assume that 



\Ei"^-DfJ/D^^J\^oo (49) 

for J7 I— > oo. This implies that the imaginary part of the 
interpolating self-energy 

--ImJ^daiE + lO) = X 

TT TT\ba\ 

-lfe.lIniSif^^(i^) ^ 

(l - 5.Re4f ^^S))' + (-|5.|Im4f ^HS) 

(50) 

has a 6 singularity at i? = E^"^ . Using Eq. (^) we can 
derive the asymptotic position 

Ei°^ ^Bd-a-l^ + Uil-{nd^a)) (51) 
and the asymptotic weight of the 5 function: 

^-lmJ:daiE + tO) ^ U^nd-a){l-{nd-a))S{E~Ej°^). 

TT 

(52) 

The weight turns out to be equal to the full weight C^^-* 
of —luiT,da{E)/Tr. Subject to the assumption (E^), Eq. 



( [5^ ) tells us that the weight of the S peak at E'^"-' will 
dominate Im'ScicriE) eventually. The real part can be 
obtained from a Kramers-Kronig-type relation: 



- -V 

TT 



lmJ:d„{E' + iO) 

E-E' 



dE' . 
(53) 

f 2") 

It turns out that it is given by Y^\J{E), which is just the 
self-energy of the SDA [cf. Eq. @]. 

Eq. ( pO[ ) shows that the imaginary part of the interpo- 
lating self-energy is non- vanishing at E = E^j'^ as well as 
at those energies E where ImS^^°^^(£;) ^ 0, i. e. within 
a certain energy interval around D^^'' / d'"^^ . We can thus 
consider Yjda{E) to consist of two additive parts. The 
dominating part of the self-energy has been identified as 
being equal to the self-energy of the SDA. The remain- 
ing part does not vanish as J7 i-^ oo as can be seen when 
expanding the interpolating self-energy ( psf ) in powers 
of 1/U directly. However, compared with the SDA part, 
its weight is smaller by a factor U"^. In fact, for strong 
U the remaining part can be neglected completely, pro- 
vided that its energetic position, which is well approxi- 



mated by the center of gravity /^da ^. 
well apart from and thereby insignificant with respect 
to the states that form the lower Hubbard band (again 
we assume that {uda) = (nd-a) < 0.5, for simplicity). 

It can be concluded that for [/ i-^ oo the interpolat- 
ing self-energy T,da{E) reduces to the SDA self-energy 

(2) 

^da i-^)' if conditions are fulfilled: The first is given 
by digj) and the second one reads: 



.(SOC) 



(54) 



for U OO. Inserting Eq. (|5l| ) and the coefficients from 
Eq. (H^) into (^9[), the first condition can be rewritten: 

f/^(n,-.)(«F)(l-(n,_.>(«F))|6.|^ 



\Bd 



B 



(HF) 
d—a 



+ ]l, + Uil-2{nd^,))\ 



■ OO . 
(55) 



As L/ I— > cx) the correlation functions Bda and B^^^^ 
stay finite. Restricting ourselves to the case (rider) = 
(nd-a) < 0.5, we have fj, i—> const.. Thus we can write 
the two conditions in the form: 



\Jl„ + U{l-2{nd-a))\ 

IJla -U (Ud-a) I OO . 



(56) 

(57) 



In the following let us discuss the implications for the 
three choices that are considered for the determination 
of the parameter /I^ according to Eqs. (^-(p5|). We 
start with the case Jl^ — /i. Obviously, both conditions 
(56) and ( |57| ) are fulfilled. Therefore, the first choice 
yields the SDA self-energy in the limit U oo. 



The second choice. 



{ndo 



v(HF) ^ 



^nda 



implies /io 



U{nd-cr)- While this fulfills the first condition, it is at 
variance with the second one. As J7 i— > oo the overall en- 
ergy dependence of the self-energy is given by the SDA; 
there are, however, non-negligible modifications for en- 
ergies around E = /D^J'^, i. e. within the lower Hub- 
bard band. Compared with the SDA, the upper Hub- 
bard band is completely unaffected, especially what con- 
cerns its spectral weight and its center of gravity. This 
implies that also the weight and the center of gravity 
of the lower Hubbard band agree with the predictions 
of the SDA, since the (zeroth and the) first moment of 
the total spectral density is reproduced exactly. In con- 
trast to the SDA, however, the non-zero imaginary part 
of the self-energy leads to quasi-particle damping within 
the lower band. Via the Kramers-Kronig-type relation 
( [53| ) the quasi-particle energies will be modified, too. 

The determination of Jlcr according to the Friedel sum 
rule (or equivalently the Luttinger theorem) is more im- 
plicit. Nevertheless, the following indirect argument can 
be given: 

We take T — 0; furthermore we again restrict our- 
selves to the case {uda) — {nd-a) < 0.5. Let us first 
mention that for the general proof ||2^ of the Friedel 
sum rule for the SIAM one has to resort to various iden- 
tities that apply to Fermi liquids. In particular, one 



needs |26 



Im T,da{E + iO) - E^ for E^O. 



(58) 



Secondly, we show that the SDA is at variance with 
the sum rule ||5||. For (uda) < 0.5 the pole of the SDA 
self-energy at E = Bd-a — M — U{1 — {ud-cr}) > lies 

(2) 

outside the range of integration. Thus T,j^J{E) is real 
and 



dE 



-U^{nd-.){l-{nd-a)) 
iE + fi-Bd-a-U{l- {ud^^W 



< 
(59) 

for all energies — oo < E <0. Since ImGdaiE -I- lO) < 
for all E and lmGdcr{E + iO) < for a certain energy 
range within — oo < i? < 0, it follows that 



Im 



iO 



Gda{E) 



oo+iO 



dE 



dE>0, 



(60) 



which according to Ref. implies that the Friedel sum 
rule is not obeyed. 

Thirdly, since is fixed by imposing the validity of 
the Friedel sum rule and since the SDA implies (|60|), 
we can conclude that the interpolating self-energy must 
be different from the SDA self-energy in all cases, espe- 
cially for f/ I— > oo. Consequently, one of the two con- 
ditions ( p6| ) or (ll^) must be violated. If it is assumed 
that the first one holds, it follows that {rid-a-)^^'^^ = 
for U larger than a certain critical value. Analogously 
to the above discussion of the case = this would 
imply ImS(;o-(0) 7^ 0, a consequence that is not com- 
patible with the Friedel sum rule. Hence, we must 
have (nd-cr)^^^^ ^ implying that for J7 1— > 00 the 
second condition does not hold. Similar to the case 
{ndtr)'^^^^ = (n-dcr), it can thus be concluded that apart 
from quasi-particle damping the overall shape of the 
spectral density follows the predictions of the SDA, in 
particular, what concerns the energetic positions and the 
spectral weights of both Hubbard bands. 

Summing up, it has turned out that all three choices 
to determine Jl^ more or less make contact with a stan- 
dard strong-coupling approach (SDA) for [/ i— > 00. This 
fact provides additional justification for the interpolat- 
ing self-energy since the SDA is known to yield rather 
satisfactory results, at least on the qualitative level. 
Therefore, we believe that the theory is able to yield re- 
liable results well beyond the weak-coupling regime. Let 
us also mention that within this context it is important 
that the theory correctly accounts for the m = 3 mo- 
ment. Otherwise, we would have ended up for U '—^ 00 
with the "Hubbard-I" self-energy ^^d^{E) only. 

VII. RESULTS 

We have evaluated the theory numerically. The pro- 
cedure is described briefly in Ref. [^. The additional 
computational effort due to the inclusion of the higher- 
order correlation functions via Bda- is almost negligible, 
and thus the algorithm remains comparatively fast. The 
results being discussed in the following have been ob- 
tained for the Bethe lattice with infinite coordination 
number. The semi-elliptic Bloch-density of states has a 
finite width W. All energies are given in units of W. 
Furthermore, we choose Cd = 0. 

Within the SDA the correlation functions Bda lead to 
an additional energetic shift of the lower and the up- 
per Hubbard band. For strong U the effective shift of 
the lower —a Hubbard band is given by {n„) B„ [ p3[ . 
Fig. 1 shows the dependence of this "band shift" on the 
occupation number n for U = 2 as obtained from our 
modified IPT. Results for the three different conditions 



(33)-(|35|) are shown. In all cases we find a non-zero but 
small band shift. Except for the case 'jla = /x, the curves 
closely resemble the corresponding results of the SDA 
HI] for a d = 3 bcc lattice. This regards the absolute 
magnitude as well as the overall dependence on n. The 
calculation for /ig. = /i yields a non-zero value for Be as 
n I— > 1, while it results in B^ — for n = 1. Contrary, 
the second and the third condition predict a continuous 
dependence on n at half-filling. Another difiiculty is ob- 
served for the case fi = fj-o + 5](t(0). Below n = 0.23 
a self-consistent solution could not be found. This also 
holds true (at a slightly different n) if we set Be = 0. 

The effect of B^ on the spectral density can be seen 
in Fig. 2 where we compare the result for the complete 
theory with the result for B^ = (the approach of Ref. 
p|]). In both cases we have chosen the third condition 
( |35[ ) to determine ]!„ . The occupation number has been 
fixed at n = 0.6 where according to Fig. 1 the band 
shift is at its maximum. Qualitatively similar to the 
spectra expected for the SIAM, both curves clearly show 
up three features: the lower and the upper Hubbard 
band at « —0.5 and E ^ 2.2 and a peak around E — 
0, which is a reminiscence of the Kondo resonance being 
strongly broadened away from half-filling. We notice 
that the difference between both spectra is rather small. 
While the low-energy features are completely unaffected, 
the upper Hubbard band slightly shifts to lower energies 
when taking into account the m = 3 moment. This is 
contrary to the SDA which predicts an energetic shift 
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Fig. 1. "Band shift" {nc,)Bcr as a function of filling n = 
(n|) + (rtj) for (7 = 2 (energies are given in units of W; W: 
width of the semi-elliptic Bloch-density of states). Dotted 
line: calculation for = fJ, [first condition (^s])]. Dashed 
line: calculation with ^o- being determined by (no-)'-^^' = 
(ria) [second condition (^)]. Solid line: /i = /lo 4- 2^(0) 
[third condition (Bq)]. 



of the upper Hubbard band (with respect to the Fermi 
energy) to higher energies by an amount (1 — n)Bu > 
0. The effect can be traced back to the (imphcit) U 
dependence of the hybridization function. 

It has not yet been finaUy clarified what is the op- 
timum choice to determine the Active chemical poten- 
tial /io-. For this purpose we compare with results from 
the exact diagonalization method of Caffarel and Krauth 
1 13, We take the data from Kajueter and Kotliar 
1 23 1 for 8 sites, U = 2 and n = 0.86. Because of the fi- 
nite number of orbitals considered in the calculation, the 
resulting spectral density is not smooth. Rather than 
comparing the spectral densities directly, a comparison 
of the integrated spectral weight is more appropriate. 
This is shown in Fig. 3. The agreement of the ED result 
with the modified IPT calculation is well provided that 
condition (|3^) or ( ]35| ) is used (curves b) and c)). In both 
cases the residual discrepancies can be attributed to the 
finite system size in the ED calculation. Significant dif- 
ferences, however, are observed between the ED result 
and the calculation for Jl^r = fi (curve a)). In particular, 
spectral weight is missing around E — —0.5. 

Fig. 4 shows the corresponding spectral densities. We 
notice that there are only minor differences between the 
results for (71^)^^^^ = (n„) and fi = fj,o + 2^(0). Apart 
from the lower and the upper Hubbard band the spec- 
tra exhibit a sharp (Kondo) resonance at £^ = 0. On 
the other hand, the spectral density that is calculated 
for pt(j = /i looks completely different. One can no longer 
distinguish unambiguously between the lower Hubbard 
band and the resonance. At -E = —0.6 a minimum can 
be found. The corresponding nearly constant trend of 




Fig. 3. Integrated spectral weight as a function of energy 
for U — 2 and n = 0.86. Solid lines: result for (a) /j.^ — fi. 
(b) {n^)^™^ = (n^>, (c) (1 = ^10 + 2^(0), Dashed line (a- 
c): exact diagonalization (8 sites), from Ref. [^, slightly 
smoothed. (The vertical scale applies to (a). The curves (b) 
and (c) have been shifted constantly.) 



the integrated weight in Fig. 3 at the same energy, how- 
ever, is at variance with the ED result which predicts a 
steep increase. Furthermore, the upper Hubbard band 
is significantly shifted to higher energies compared with 
the results for the second and the third condition which 
according to Fig. 3 reliably reproduce the peak position. 

There is one unsatisfactory property of the modified 
IPT when using the second condition 
namely the drop of the spectral density just below the 
Fermi edge = as can be seen in Fig. 4. This behav- 
ior, however, is only found for very strong interaction. 
Fig. 3 shows results ior U — 1 which still means strong 
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Fig. 2. Spectral density as a function of energy for (7 = 2 
and n = 0.6. Calculations assuming jj, = fio + Ect(O). Solid 
line: complete theory. Dashed line: Ba = 0. 
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correlation. Here we notice an almost perfect agreement 
between the results for the second and the third condi- 
tion. Although much closer agreement between all three 
versions is obtained generally if the interaction is re- 
duced from U = 2 to U =1, the result for Jl„ = is 
still significantly different. Remarkably, all three curves 
in Fig. 5 almost coincide at E = Q. According to Ref. 

this indicates that Luttinger's theorem is fulfilled 
approximately when using the first (^3|) or second (jsj) 
condition. 

VIII. CONCLUSIONS AND OUTLOOK 

In this paper we have presented a modification of the 
recent approach of Kajueter and Kotliar Using the 
self-consistent mapping onto the SIAM, an approximate 
analytical expression for the self-energy of the infinite- 
dimensional Hubbard model could be constructed that 
reproduces a number of exactly solvable limits. The con- 
ceptual improvement consists in the consideration of an 
additional, the m — 3, moment of the spectral density. 
It has been shown that the higher-order correlation func- 
tions that are included in the extra term can be ex- 
pressed without further approximations by means of the 
spectral density. This allows for a self-consistent (nu- 
merical) solution. The additional computational effort 
needed is almost negligible. The Green function on the 
real axis at T = can be computed fast compared with 
QMC or ED techniques. 

The theory contains a fictive chemical potential Jla- 
which is considered as a free parameter that can be 
fixed by a rather arbitrary condition without losing rigor 
in all limiting cases mentioned. In this paper we have 
taken into account three different possibilities to deter- 
mine Jlcr- The numerical results prove that the different 



choices may imply considerable differences between the 
shapes of the resulting spectral densities, especially for 
very strong interaction U. Thus further information was 
needed to get a conclusive theory. We have compared 
our results with the data of an ED study taken from Ref. 

. By the comparison the most simple choice Jla = fi 
is excluded. 

For both choices, ^ = + ^£,(0) and (ricr)^^^^ — 
(ricr), we encountered a minor difficulty: self-consistent 
solutions could not be found for fillings below n = 0.23 
(at U — 2) when taking the first one; using the latter, 
we observed an implausible drop of the spectral density 
just below the Fermi edge which, however, is present in 
the case of very strong interaction only. 

The condition /z = yLto-l-So-(O) represents the Luttinger 
theorem for the d = oo Hubbard model. Imposing the 
Luttinger theorem as a condition to fix Jla- as has been 
suggested by Kajueter and Kotliar |^ implies a consid- 
erable restriction of the theory: in this form the theo- 
rem is only meaningful for a paramagnet at T = [ pO[ . 
This disadvantage is not present when using the condi- 
tion which was introduced by Martin- 
Rodero et al. originally ^ . Finite temperatures and 
ferro- or antiferromagnetism can be treated without dif- 
ficulties. Furthermore, the condition is much easier to 
handle numerically. On the other hand, the difference 
found between the numerical results for the spectral den- 
sity using either ^ = fj,Q -\- I]o-(0) or (no-)^^^^ = (ria) are 
rather small, and the agreement with the ED data is 
equally good. 

The usefulness of the m = 3 moment is apparent in 
the limit of strong correlations U t-^ oo. This limit of 
the approach has been investigated within the SIAM, i. 
e. for a fixed hybridization function A(i?). The mean 
energetic positions and the weights of the upper and the 
lower Hubbard band agree with the predictions of the 
SDA and with the exact results of Harris and Lange 
p^ . Here the m = 3 moment turns out to be decisive. 
Otherwise, one would have ended up with the "Hubbard- 
I" solution only. 

The results for the paramagnetic d ^ oo Hubbard 
model on the Bethe lattice at T = have shown the ef- 
fect of Ba on the spectral density to be rather small. 
Previous studies, however, strongly suggest that the 
m = 3 moment is quite important in the context of spon- 
taneous magnetism. This is obvious, for instance, when 
comparing the SDA (correct moments up to to = 3) 
with the Hubbard-I solution (correct moments up to 
m = 2). While the Hubbard-I solution yields magnetic 
order only under extreme circumstances, magnetism is 
favored within the SDA: The term B^ opens the possibil- 



ity for a spin-dependent band shift. Consistent with the 
results found here, the effect of B^- in the paramagnetic 
phase is small within the SDA as well |5^. Comparing 
the Hubbard-III alloy-analogy solution with a recently 
developed modification |^ where again B^r is included 
additionally, also stresses the importance of the to = 3 
moment for spontaneous magnetism. 

The application of the presented method to magnetic 
phases represents an interesting task for future studies. 
Let us mention that ferromagnetism in the d = oo Hub- 
bard model for an fcc-type lattice has been found re- 
cently in a QMC calculation Particle-hole sym- 
metry requires = at half-filling (for a symmetric 
Bloch-density of states) in the paramagnetic phase. In 
this case the usual IPX is recovered. However, B^ ^ 
is possible for an antiferromagnet at half-filling. Future 
work may thus check whether the approach can improve 
the IPX results for antiferromagnetic order at rt = 1 
which are not completely satisfactory pi] ]. 
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